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ABSTRACT 

Understanding cosmic reionization requires the identification and characterization of early sources of 
hydrogen-ionizing photons. The 2012 Hubble Ultra Deep Field (UDF12) campaign has acquired the deep- 
est infrared images with the Wide Field Camera 3 aboard Hubble Space Telescope and, for the first time, 
systematically explored the galaxy population deep into the era when cosmic microwave background (CMB) 
data indicates reionization was underway. The UDF12 campaign thus provides the best constraints to date on 
the abundance, luminosity distribution, and spectral properties of early star-forming galaxies. We synthesize 
the new UDF12 results with the most recent constraints from CMB observations to infer redshift-dependent 
ultraviolet (UV) luminosity densities, reionization histories, and electron scattering optical depth evolution con- 
sistent with the available data. Under reasonable assumptions about the escape fraction of hydrogen ionizing 
photons and the intergalactic medium clumping factor, we find that to fully reionize the universe by redshift 
z ~ 6 the population of star-forming galaxies at redshifts z ~ 7 - 9 likely must extend in luminosity below the 
UDF12 limits to absolute UV magnitudes of Muv ~ -13 or fainter. Moreover, low levels of star formation 
extending to redshifts z ~ 15-25, as suggested by the normal UV colors of z ~7-8 galaxies and the smooth 
decline in abundance with redshift observed by UDF12 to z — 10, are additionally likely required to reproduce 
the optical depth to electron scattering inferred from CMB observations. 
Subject headings: cosmology:reionization - galaxies:evolution - galaxies:formation 



1. INTRODUCTION 

The process of cosmic reionization remains one of the most 
important outstanding problems in galaxy for mation and cos- 
molo gy. After recombination at z ~ 1090 (Hinsha w et alj 
12012ft . gas in the univer se was mostly n eutral. How- 
ever, observations of the |Gu nn & Peterson 
in quasar spectra (e.g. , iFan et al.l 120011 120021 



1965b trough 
2001 [2006b; 

Diorg ovski et al.ll2.001h indicate that intergalactic gas has be- 
come almost fully reionized by redshift z ~ 5. The elec- 
tron scattering optical depth inferred from CMB observations 
suggests that if the universe was instantaneously reionized, 
then reionization would occur as early a s redshift z ~ 10 
( Sper gel et all 12003b iHinshaw et all 120121) . Given the dra- 
matic decline in the abundance of quasars beyond redshift 
z ~ 6, they very likely ca nnot be a significant contributor to 
cosm ic reionization (e.g.. IWillott et al.l 120 1 0b iFontanot et al.l 
120121) eve n though quasars hav e been discovered as early 
as z ~ 7 (Mortlock et al] 1201 ll) . Star-forming galaxies at 
redshifts z > 6 have therefore long been postulated as the 
likely agents of cosmic reionization, and their time-dependent 
abundance and spectral properties are thus crucial ingredients 
for understanding how intergalactic hydrogen became reion- 
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ized (for reviews, s ee jFan et al.ll2006al iRobertson et al.ll20lol 
lLoeb & Furlanettoll2012l) . 

We present an analysis of the implications of the 2012 Hub- 
ble Ultra Deep Fielcfj (UDF12) campaign results on the abun- 
dance and spectral characteristics of galaxies at z ~ 7- 12 for 
the reionization process. The UDF12 campaign is a 128-orbit 
Hubble Space Telescope (HST) program (GO 12498, PI: El- 
lis) with the infrared (IR) channel on the Wide Field Cam- 
era 3 (WFC3/IR) that acquired the deepest ever images in the 
IR with HST in the Hubble Ultra Deep Field (HUDF) in Fall 
2012 (the UDF1 2 p roject and data overview are described in 
lEllis et al.ll20L3l and lKoekemoer et al.ll2013l) . Combined with 
previous HUDF observations (GO 11563, PI: G. Illingworth; 
GO 12060, 12061, 12062, Pis: S. Faber and H. Ferguson; GO 
12099, PI: A. Riess), the UDF12 imaging reaches depths of 
F105 = 30, J 125 = 29.5, J m = 29.5, and H m = 29.5 (5 - a AB 
magnitudes). The UDF12 observations have provided the first 
det erminations of th e galaxy abundance at redshifts 8 .5 < z < 
12 (lEllis et alJl20T3l) . precise determinati ons of the galaxy lu- 
minosity function at redshifts z ~ 7-8 dSchenker et aT] 2012a; 
iMcLure et al.l 120121) . robust ultraviolet (UV) spe ctral slope 
measurements at z ~ 7-8 dDunlop et al.N2012bl) , and size- 
lumino sity relation measures at redshifts z ~ 6- 8 dOno et alj 
l2012ah . 

Our earlier UDF12 publications already provide some new 
constraints on the role that galaxies play i n cosmic r e ioniza - 
tion and the duration of the process. In Ellisetal.1 d20ll . 
we argued that continuity in the declining abundance of star- 
forming galaxies over 6 < z < 10 (and possibly to z — 12), 
implied the likelihood of further star formation beyond the 
redshift limits currently probed. Likewise, in iDunlop et alj 
(2012b), the constancy of the UV continuum slope measured 
in z — 7-9 galaxies over a wide range in luminosity supports 

8 http://udfl2.arizona.edu 



2 



Robertson et al. 



the contention that the bulk of the stars at this epoch are al- 
ready enriched by earlier generations. Collectively, these two 
results support an extended reionization process. 

We synthesize these UDF12 findings with the recent 9- 
year Wilkinson Micro wave Anisotropy Probe (WMAP) results 
(Hinshaw et al. 2012) and stellar mass density measurements 
dStark et al.ll2012l) to provide new constraints on the role of 
high-redshift star-forming galaxies in the reionization pro- 
cess. Enabled by the new observational findings, we perform 
Bayesian inference using a simple parameterized model for 
the evolving UV luminosity density to find reionization his- 
tories, stellar mass density evolutions, and electron scatter 
optical depth progressions consistent with the available data. 
We limit the purview of this paper to empirical modeling of 
the reionization process, and comparisons with more detailed 
galaxy formation models will be presented in a companion 
paper (Dayal et al., in preparation). 

Throughout this paper, we assume the 9-year WMAP cos- 
mological parameters (as additionally constrained by exter- 
nal CMB datasets; h = 0.705, fi m = 0.272, fi A = 0.728, 
fib = 0.04885). Magnitu des are reported using the AB sys- 
tem dOke & Gunnlll983l) . All Bayesian inference and maxi- 
mum likelihood fitting is performed usin g the MultiNest code 
dFeroz & Hobsoril2008t lFeroz et al.ll2009b . 

2. THE PROCESS OF COSMIC REIONIZATION 

Theoretical models of the reionization process have a long 
history. Early analy tic and numerical models of the reioniza- 
tion process (e.g.. IMadau et al.l [l999t [Miralda-Escu de et al. 
2000[ lGnedinll2000t iBarkana & LoebfcoqU iRazoumov et al. 
2002t IWvTthe & Loebl 120031: iCiardi et all 120031) highlighted 
the essential physics that give rise to the ionized inter- 
galactic medium (IGM) at late times. In the follow- 
ing description of the cosmic reioniza tion process we fol - 
low most closely the mo de ling of IMadau et alj d 19991) . 
Bolton & Haehneltl d2007bl) iRobertson et all d2010l) , and 
Kuhlen & Faucher-Giguere (2012), but there has been closely 
related recent work by Ciardi et al.l (12012b and Uensen et al.l 
d20Tl . 

The reionization process is a balance between the recombi- 
nation of free electrons with protons to form neutral hydrogen 
and the ionization of hydrogen atoms by cosmic Lyman con- 
tinuum photons with energies E > 13.6 eV. The dimension- 
less volume filling fraction of ionized hydrogen Qmi can be 
expressed as a time-dependent differential equation capturing 
these competing effects as 



where dotted quantities are time derivatives. 
The comoving density of hydrogen atoms 

{nn)=X v Q, b p c (2) 

depends o n the primordia l mass-fraction of hydrogen X p = 
0.75 (e.g.. lHouetaT1l2011l) . the critical density p c = 1.8787 x 
10~ 29 /T 2 g cm 3 , and the fractional baryon density fib- 
As a function of redshift, the average recombination time in 
the IGM is 

W = [C mi a B (T)(l +F p /4Z p )(n H )(l +z) 3 ]"' . ( 3 ) 

where a^(T) is the case B recombination coefficient for hy- 
drogen (we assume an IGM temperature of T = 20,000K), 
Y p = 1 -X p is the primordial helium abundance (and accounts 



for the numbe r of free electrons per proton in th e fully ion- 
ized IGM, e.g., Kuhlen & Faucher-Giguere 2012), and Chii = 
(n H ) / («h) 2 is the "clumping factor" that accounts for the ef- 
fects of IGM inhomogeneity through the quadratic depen- 
dence of the recombination rate on density. 

Simulations suggest that the clumping fact or of IGM gas 
is Chj t ~ 1 — 6 at the redshifts of interest (e . g.. ISokasian et al.l 
poTlIliev et al.M2006h iPawlik et al.ll2009t iShull et al.li2012h 
iFinlator et al.ll2012l) . While early hydrodynamical simulation 
studies suggested that the clumping fact or could be as high 
as Ch tt ~ 10-40 at redshifts z < 8 (e.g.. iGnedin & Ostrikerl 
1997), recent studies that separately identify IGM and inter- 
stellar medium gas in the simulations and employ a more de- 
tailed modeling of the evolving UV background have found 
lower values of Chii- An interesting study of the redshift evo- 
lution of the clumping factor was provided by Pawli k et al.l 
(2009, see their Figures 5 and 7). At early times in their 
simulations (z > 10), the clumping factor was low (Chii < 3) 
but increased wit h decreasing redshift at a rate similar to pre- 
dictions from the Miralda-Esc ude et"ai1 (120001) model of the 
evolving IGM density probability distribution function. In 
the absence of photoheating, at lower redshifts (z < 9) the 
clu mping factor would b e gin to increase more rapidly than 
the Miralda-Es cude et al.l (120001) prediction to reach Chii ~ 
10-20 by z ~ 6. In the presence of photoheating, the evo- 
lution of the clumping factor depends on the epoch when the 
uniform UV background becomes established. If the IGM 
was reheated early (z ~ 10-20) the predicted rise in clump- 
ing factor breaks after the UV background is established and 
increases only slowly to Che ~ 3 — 6 at late times (z ~ 6). If 
instead, and perhaps more likely, the IGM is reheated later 
(e.g, z ~ 6-8) the clumping factor may actually decrease at 
late times from Chii ^6-10atz~8to C hii ~ 3-6 at z ~ 6. 

The results of th ese simulations (e.g.. IPawlik et al.l 120091: 
Finlatoret al] 120121) in part motivate our choice to treat the 
clumping factor as a constant Chii ~ 3 since over a wide range 
of possible redshifts for the establishment of the UV back- 
ground the clumping factor is expected to be C hii ^2-4 
at z < 12 (see Figure 5 of Paw lik et al.l 12009) and lower 
at earlier times. In comparison with our previous work 
dRobertson et al.ll2010l) . where we considered Chii = 2-6 and 
frequently used Chii = 2 in Equation [3] we will see that our 
models complete reionization somewhat later when a some- 
what larger value of Chii is more appropriate. However, 
we note that the end of the reionization process may be 
more complicated th an what we have descri bed above (see, 
e.g., Section 9.2.1 of lLoeb & Furlanettoll2012l) . As reioniza- 
tion progresses, the ionized phase penetrates more and more 
deeply into dense clumps within the IGM - the material that 
will later form the Lyman-a forest (and higher column den- 
sity systems). These high-density clumps recombine much 
faster than average, so Chi i may increase throughout reioniza- 
tion (|Furlanetto & Ora 120051) . Combined with the failure of 
EquationQ~|to model the detailed distribution of gas densities 
in the IGM, we expect our admittedly crude approach to fail 
at the tail end of reionization. Fortunately, we are primarily 
concerned with the middle phases of reionization here, so any 
unphysical behavior when gnu is large is not important for 
us. 

The comoving production rate n; on of hydrogen-ionizing 
photons available to reionize the IGM depends on the intrin- 
sic productivity of Lyman continuum radiation by stellar pop- 
ulations within galaxies parameterized in terms of the rate of 



UDF12: New Constraints on Cosmic Reionization 



3 



hydrogen-ionizing photons per unit UV (1500A) luminosity 
£km (with units of ergs" 1 Hz), the fraction / esc of such photons 
that escape to affect the IGM, and the total UV luminosity 
density puv (with units of ergs s" 1 Hz -1 Mpc~ 3 ) supplied by 
star-forming galaxies to some limiting absolute UV magni- 
tude Muv- The product 

"ion = /esc6onPUV (4) 

then determines the newly available number density of Ly- 
man continuum photons per second capable of reionizing in- 
tergalactic hydrogen. We note that the expression of n\ oa in 
terms of UV luminosity de nsity rather than star formation rate 
(c.f jRobertson et aLll2010t) is largely a matter of choice; stel- 
lar population synthesis models with assumed star formation 
histories are required to estimate £j on and using the star for- 
mation rate density ps FR in Equation [4] therefore requires no 
additional assumptions. Throughout this paper, we choose 
/e S c = 0.2. As shown by Ouch i et aTl ([2009), escape fractions 
comparable to or larger than / esc = 0.2 during the reionization 
epoch are required for galaxies with typical stellar populations 
to contribute significantly. We also consider an e volvi ng / esc 
with redshift, with the results discussed in Section l6T2l below. 

The advances presented in this paper come primarily from 
the new UDF12 constraints on the abundance of star-forming 
galaxies over 6.5 < z < 12, the luminosity functions down to 
Muv — —17, and robust determinations of their UV continuum 
colors. For the latter, in Section |3] we use the UV sp ectral 
slope of high-redshift galaxies by iDunlop et al.l (f2 012b) and 
the ste llar population synthesis models of Bruzual & Chariot 
( 120031) to inform a choice for the number £j on of ionizing pho- 
tons produced per unit luminosity. For the former, the abun- 
dance and lu minosity di s tributi o n of high-redshi f t galax ies de- 
termined by lEllis et al.l d2013l) . ISchenker et al.l (l2012al) . and 
iMcLure et al.1 (120121) provide estimates of the evolving UV 
luminosity density puv- The evolving UV luminosity density 
supplied by star-forming galaxies brighter than some limiting 
magnitude Muv is simply related to an integral of the lumi- 
nosity function as 

/Muv 
<5>{M)L(M)dM, (5) 
oo 

where L is the luminosity and the functional fo rm of the 
galaxy luminosity function is often assumed to be a lSchechter] 
(119761) function 

$(M) = 0.4 In 10 fa [lO a4(M *- M) ] 1+Q exp [-10°- 4 ( M *" M >] (6) 

parameterized in terms of the normalization fa (in units 
of Mpc" 3 mag" 1 ), the characteristic galaxy magnitude M t , 
and the faint-end slope a. Each of these parameters may 
evolve with redshift z, which can affect t he relative impor- 
tance of faint galaxies f or reionization (e.g., O esch et al.l2009t 
iBouwens et al.l l2012bl) . In Sections [4] and 14.11 below, we 
present our method of using the previous and UDF12 data sets 
to infer constraints on the luminosity density as a function of 
redshift and limiting magnitude. 

2.1. Stellar Mass Density as a Constraint on Reionization 

The UV luminosity density is supplied by short-lived, mas- 
sive stars and therefore reflects the time rate of change of the 
stellar mass density p+(z). In the context of our model, there 
are two routes for estimating the stellar mass density. First, 
we can integrate the stellar mass density supplied by the star 



formation rate inferred from the evolving UV luminosity den- 
sity as 

pjz) = (1 -R)J~ ihfAz') P vv(z')^dz\ (7) 

where ry s f r (z) provides the stellar population model-dependent 
conversion between UV luminosity and star formation rate (in 
units Mq yr" 1 ergs" 1 s Hz), R is the fraction of mass returned 
from a stellar population to the ISM (28% for a lSalpeteril 19551 
model after ~ 10 Gyr for a 0.1 - 100M© IMF), and dt/dz 
gives the rate of change of universal time per unit redshift 
(in units of yr). While rj^iz) is in principle time-dependent, 
there is no firm evidence yet of its evolution and we adopt a 
constant value throughout. 

While the evolving stellar mass density can be calculated 
in our model, the observational constraints on have in- 
volved integrating a composite stellar mass function con- 
structed from the UV lumi nosity function and a stellar mass 
to UV luminos ity relation (jGonzalez et al.l 1 201 U IStark et aT] 
l20H see also lLabbe et al.l 120121) . Using near-IR observa- 
tions with the Spitzer Space Telescope, stellar masses of UV- 
selected ga l axies are measured as a function of luminosity. 
IStark et aU (120121) find that the stellar mass - Luv relation, 
corrected for nebular emission contamination, is well de- 
scribed by 

logm* = 1.4331ogL uv -31.99 + .A(z) (8) 

where m* is the stellar mass in Mq, Luv is the UV luminosity 
density in ergs s" 1 Hz" 1 , and f+(z) = [0,-0.03,-0.18,-0.40] 
for redshift z« [4,5,6,7]. 

The stellar mass density will then involve an integral over 
the product of the UV luminosity function and the stellar mass 
m^(Luv)- However, the sign ificant scatter in the m^(Luy) rela- 
tion ( a ft 0.5 in logm*, see Gonz alez et alJl2QTlT : IStark et al.l 
120121) must be taken into account. The gaussian scatter p[M'— 
M(logm+)] in luminosity contributing at a given logm* can be 
incorporated into the stellar mass function dn+/dlogm+ with 
a convolution over the UV luminosity function. We write the 
stellar mass function as 

dn * = dM [ <S>{M')p[M'-M{\ogm+)]dM'. (9) 
dlogm* dlogm+J^ 

For vanishing scatter Equation [9] would give simply 
dn+jdm* = $[M(m*)] x dM/dm+. The stellar mass density 
p± can be computed by integrating this mass function as 

f logm * dn* 

pj<m in z) = / — -m'JXogm'*. (10) 

A primary feature of the stellar mass function is that the stellar 
mass - UV luminosity relation and scatter flattens it relative 
to the UV luminosity function. Correspondingly, the stellar 
mass density converges faster with decreasing stellar mass or 
luminosity than does the UV luminosity density (Equation|5]). 
The stellar mass density then serves as an additional, integral 
constraint on ri\ on . 

2.2. Electron Scattering Optical Depth 

Once the evolving production rate of ionizing photons h mn 
is determined, the reionization history Qmiiz) of the universe 
can be calculated by integrating Equation [TJ An important 
integral constraint on the reionization history is the electron 
scattering optical depth r inferred from observations of the 
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CMB. The optical depth can be calculated from the reioniza- 
tion history as a function of redshift z as 

r(z)= rc(n H )<j T f e Q ml (z')H(z')(l+z') 2 dz\ (11) 
Jo 

where c is the speed of light, ot is the Thomson cross sec- 
tion, and H{z) is the redshift-dependent Hubble parameter. 
The number / e of free electrons per hydrogen nucleus in 
the ionize d IGM depends on the ionization state of helium. 
Following Kuh len & Faucher-Giguerel (12.012b and other ear- 
lier works, we assume that helium is doubly ionized (f e = 
1 +Yp/2X P ) at z < 4 and singly ionized (f e = 1 + Y P /4X P ) at 
higher redshifts. To utilize the observational constraints on 
r as a constraint on the reionization history, we employ the 
posterior probability distribution p(f) determined from the 
Monte Carlo Markov Chains used in the 9-year WMAP re- 
sultfl as a marginalized likelihood for our derived r values. 
This method is described in more detail in Section [5] 

3. UV SPECTRAL SLOPES AND THE IONIZING PHOTON BUDGET 

A critical ingredient for determining the comoving produc- 
tion rate of hydrogen ionizing photons is the ratio ^ on of the 
Lyman continuum photon emission rate per unit UV (1500A) 
luminosity spectral density of individual sources. Since Ly- 
man continuum photons are predominately produced by hot, 
massive, UV-bright stars, it is sensible to expect that £j on will 
be connected with the UV spectral slope of a stellar popula- 
tion. Here, we use observational constraints on the UV slope 
of high- redshift galaxies deter mined from the UDF12 cam- 
paign by Dunlo pet alj d2012bl) and stell ar population synthe- 
sis models by Bruzual & Chariot (2003, BC03) to estimate a 
physically-motivated £j on consistent with the data. In this pa- 
per, we concentrate on placing constraints on £; on ; for a much 
more detailed analysis and int erpretation of the UV slope re- 
sults from UDF12, please see Dunlop et al. C26 l2bb . 

Prior to the UDF12 program, observations of high-redshift 
galaxies in the HUDF09 WFC3/IR campaign provided a first 
estimate of the UV spectral slopes (/?, where f\ oc A") of 
z > 7 galaxies. Early results from the HUDF09 team in- 
dicated that these high-redshif t galaxies had extrao rdinar- 
ily blue UV slopes of /3 w — 3 dBouwens et al.ll2010h . much 
bluer than well-studied st arburst galaxies at lower redshifts 
(e.g., iMeurer et al.l Q999). In the intervening period be- 
fore the UDF12 data was acq uired, several workers argued 
against such e xtrem e values (Finkelste in et al.l 120101 120121: 
McLu reet al.l [2QlTt IDunlop et al.l l2012at iBouwens et alJ 
2012at [Rogers et alj|2012l) . The UDF12 campaign provided 
significantly deeper H\(,o imaging data used in the spectral 
slope determination (e.g., j3 = 4.43(7i25 - Hi6o) ~~ 2) at red- 
shifts z ~ 7 - 8) and added /140 imaging that reduces poten- 
tial observational biases and enables a first UV slope deter- 
mination at z ~ 9. T hese measurements were presented by 
IDunlop et al.l d2012bl) . whose results are discussed in the con- 
text of the present pa per in Figure [TJ (data points, left panel). 
IDunlop et al.l d2012b) measured the spectral slope /3 as a func- 
tion of galaxy luminosity and redshift in the range -19.5 < 
Muv <-17.5atz^7-8. Using their reported measurements 
(see their Table 1), we performed simple fits of a constant to 
their j3 values at each redshift separately and found maximum 
likelihood values of /3(z ~ 7) = -1.915 and /3(z - 8) = -1 .970 
(Figure [TJ red lines in left panel) consistent with the single 

9 http://lambda.gsfc.nasa.gov 



M uv =-18 z ~ 9 measurement of /3(z - 9) = -1.80 ±0.63. 
The 68% credibility intervals on a constant /3 at each redshift 
(Figure [TJ grey areas in left panel) suggest that across red- 
shifts z ~ 7 - 9 galaxies are consistent with a non-evolving 
UV spectral slope in the range -2.1 < j3 < —1.7. The appar- 
ent constancy of j3 with redshift avoids the need for strong 
assumptions about the redshift evolution of galaxy properties. 

To connect these UV spectral slope determinations to a 
value of £i on , we must rely on stellar population synthesis 
models. We use the standard BC03 models to extract model 
spectral of stellar populations with a range of star formation 
histories (bursts and constant star formation rates), metallic- 
ities (Z = 0.0001 -0.05), dust absorptio n (A v w 0.1- 1; cal- 
culated using the Chariot & Fall 2000 dust model) and ini - 
tial mass functions (IMF: IChabrierfl2003l and lSalpeterill955l) . 
These models provide both the Lyman continuum produced 
by hot stars and the full spectral energy distributions (SEDs) 
of the composite stellar population per unit star formation rate 
(SFR) or stellar mass. We determine for each model by di- 
viding the Lyman continuum photon production rate per unit 
SFR or stellar mass by a 1500A luminosity spectral density 
(in ergs s" 1 Hz -1 ) per unit SFR or stellar mass measured with 
a synthetic filter with a flat response and 100A width. The 
model SFR or stellar mass then scales out of the ratio £i on 
(with units ergs" 1 Hz). By measuring the UV slope values 
of each model using synthetic photometry with the J\25 and 
Hi(,o total throughput response on SEDs redshifted appropri- 
ately for z ~ 7 observations, the value of £; on for each model 
j3 can be studied as a function of age, metallicity, IMF, and 
star formation history. We have checked that our methods for 
measuring j3 and IR colors from the synth etic model spec- 
tra re produce results from the literature (see Robertso n et alJ 
120071 and Figure 2 of lRogers et aDl2012l) . 

The right panel of Figure [TJ shows the £i on of a variety of 
BC03 constant SFR models as a function of UV slope 0, 
over the range -2. 1 </3 <— 1.7 suggested by the lDunlop et aTJ 
d2012bl) color measurements and further constrained such that 
the age of the stellar populations is less than the age of the 
universe at redshift z ~ 7 (t sa 7.8 x 10 s yr). We will limit our 
further discussion of BC03 models to constant SFR histories, 
as we find single bursts populations display too wide a range 
of £km with j3 to be tightly constrained by /3 alone (although 
limits on the luminosity will constrain the available £j on for a 
given j3 for single burst models, without the individual SED 
fits to objects we cannot usefully constrain such models over 
the wide range of luminosity and redshifts we examine). 

Three broad types of BC03 constant SFR models are con- 
sistent with values of j3 = -2. Genetically, the BC03 con- 
stant SFR models evolve from large values of £j on and large 
negative values of j3 at early times to a roughly horizontal 
evolution with constant £j on with f3 increasing at late times 
(t > 10 8 yr). Metal poor (Z < Z Q ) constant SFR populations 
without dust produce values of j3 sig nificantly bluer (mor e 
negative) than observed at z ~ 7-9 (Dunlop et al. 2012b). 
Mature (> 10 s yr old), metal-rich (Z ~ Z ), dust free stel- 
lar populations evolve to a constant value of log£j on ps 24.95 - 
25.2 log ergs -1 Hz over the observed j3 values (as Figure [TJ 
indicates, the results are largely independent of IMF for con- 
stant SFR models). Applying a dus t absor ption of Ay ~ 0.1 
using the model of Chariot & Falll d2000l with parameters 
Ty = 0.25 and an ISM attenuation fraction of 0.3) shifts the 
SED evolution tracks down in £j on (from dust absorption) and 
to redder /3, such that young, metal-rich stellar populations 
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FIG. 1 . — Spectral properties of high-redshift galaxies and the corresponding properties of stellar populations. Dunlop et al. (2012b) used the new UDF12 HST 
observations to measure the UV spectral slope ft of z ~ 7— 9 galaxies as a function of luminosity (data points, left panel). As the data are consistent with a constant 
ft independent of luminosity, we have fit constant values of ft at redshifts z ~ 7-8 (maximum likelihood values of ft(z ~ 7) = —1.915 and ft(z ~ 8) = -1.970 shown 
as red lines, inner 68% credibility intervals shown as grey shaded regions; at z ~ 9 the line and shaded region reflect the best fit value of 0(z ~ 9) = — 1 .80 ± 0.63). 
The data are broadly consistent with ft = —2 (indicated with grey band in right panel), inde pendent of re dshift and luminosity. To translate the UV spectral slope 
to a ratio £i on of ionizing photon production rate to UV luminosity, we use the Bruzual & Chariot ( 2003, BC03) stellar population synthesis models (right panel) 
assuming a constant star formation rate (SFR). The constant SFR models evolve from a declining with increasing ft at early times to a relatively flat £; on at 
late times (we plot the values of £j on vs. ft for population ages less than the age of the universe at z ~ 7, t = 7.8 X 10 s yr). Three broad types of BC03 constant 
SFR models are consistent with value s of ft = —2: mature (> 10 8 yr old), metal-rich, dust free stellar populations, mature, metal-rich stellar populations with 
dust (Ay ~ 0.1 calculated using the Chariot & Fall (2000) model), and young , metal-rich stellar populations with dust. Dust free models are plotted with solid 
lines, while dusty models are shown as dashed lines. We assume the Chabrier (2003) initial mass function (IMF), but the Salpeter (1955) IMF produces similar 
values of £i on (dotted lines, dust-free case shown). Based on these models we optimistically assume log£j on = 25.2 log ergs - ' Hz, but this value is conservative 
compared with assumptions widely used in the literature. 



with dust can also reproduce values of j3 « — 2 while maintain- 
ing log£i on = 24.75-25.35 log ergs" 1 Hz. Moderately metal- 
poor models (Z ~ 0.2-0.4Z©) with as much as Ay ~ 0.1 can 
also reproduce /3 « — 2 for population ages t > 10 8 yr, but the 
most metal poor models in this range require t > 4 x 10 8 yr (an 
initial formation redshift of z > 12 if observed at z ~ 7). We 
note that our conclusions about the connection between f3 (or 
J\25~H\qq color) and the properties of stellar populations is 
wholly con sistent with previous re sults in the literature (e.g., 

Figure 7 of iFinkelstein et al.ll2010h. 

While the (3 measurements of Dunlop et al.l (12012b!) have 
greatly reduced the available BC03 stellar population model 
parameter space, there is still a broad allowable range of 
log Cion ~ 24.75-25.35 log ergs" 1 Hz available for constant 
SFR models with UV spectral slopes of (3 ~ -2. We therefore 
adopt the value 



log 6on = 25.2 log ergs Hz (adopted) 



(12) 



throughout the rest of the paper. This is in the upper 
range of the available values shown in Figure Q] but is com- 
parable to values adopted elsew here in the literature (e.g., 
iKuhlen & Faucher<j iguere 2012, who assume log w25.3 
for (3 = -2, see their Equations 5 and 6). 

We also considered stellar p opulations reddened by neb- 
ular continuum emission ("e.g.. [Schaerer & de BarrosTl2009l 
I2010t lOno et alj|20!ot iRobertson et al.ll2010l) , which in prin- 
ciple could allow relatively young, metal poor populations 
with larger £j n to fal l into the window of (3 values found by 
iDunlop etall d2012bl) . We find that for / esc ~ 0.2, nebular 



models applied to young (< 100 Myr) constant star formation 
rate BC03 models are still marginally too blue (f3 ~ -2.3). Al- 
though more detailed modeling is always possible to explore 
the impac t of nebular emissio n on £i on , the uniformity ob- 
served by Dunlop et al. (2012b) in the average value of (3 over 
a range in galaxy luminosities may argue against a diverse 
mixture young and mature ste llar populations in th e current 
z — 7-8 samples. However, as Dunlop et al. (2012b) noted, a 
larger intrinsic scatter could be present in the UV slope distri- 
bution of the observed population but not yet detected. Simi- 
larly, top heavy initial mass function stellar populations with 
low metallicity like the 1 - 100M(7) Salpeter I MF models of 
ISchaereri (120031) used bv lBouwens et al.l (I2010h to explain the 
earlier HUDF09 data, are disfavored owing to their blue spec- 
tral slopes. 

For reference, for conversion from UV luminosity spectral 
density to SFR we note that for po pulation ages t > 10 s yr a 
constant SFR BC03 model with a IChabrierl d2003b IMF and 
solar metallicity provides a 1500A luminosity spectral density 
of 



^uv 



1.25 x 10 28 x 



SFR 

M Q yr" 



ergs s Hz 



(13) 



while, as noted bv lMadau et alj d 1998b . a comparable Salpeter 
(119551) model provides 64% of this UV luminosity. A very 
metal -poor population (Z = Zq/200) would provide 40% 
more UV luminosity per unit SFR. 

4. ULTRAVIOLET LUMINOSITY DENSITY 
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In addition to constraints on t he spectral energy distri- 
butions of high-redshift galaxies (Dun lop et al.ll2012bl) . the 
UDF12 observations provide a critical determination of the 
luminosity function of star forming galaxies at redshifts 7 < 
Z < 9. As described in Section [2] when calculating the co- 
moving production rate n loa of hydrogen ionizing photons per 
unit volume (Equation|4|l the UV luminosity density puv pro- 
vided by an integral of the galaxy luminosity function is re- 
quired (Equation 01. An accurate estimate of the puv pro- 
vided by galaxies down to observed limits requires a care- 
ful analysis of star-forming gal axy samples at faint m agni- 
tudes. Using the UD F12 data, iSchenker etafl (I2012al) and 
iMcLure et al.l (|2012) have produced separate estimates of 
the z ~ 7 - 8 galaxy luminosity function for different sam- 
ple selections (color-selected drop-out and spectral energy 
distribution-fitted samples, respectively). As we demonstrate, 
the UV luminosity densities computed from these separate lu- 
minosity functions are consistent within 1 - a at z ~ 7 and in 
even closer agreement at z ~ 8. Further. iMcLure et al.l (120 1 2b 
have provided the first luminosity function estimates at z ~ 9. 
Combined, these star-forming galaxy luminosity function de- 
termination provide the required constraints on puv in the 
epoch z > 7 when, as we show below, the ionization fraction 
of the IGM is likely changing rapidly. 

Given the challenge of working at the limits of the ob- 
servational capabilities of HST and the relatively small vol- 
umes probed by the UDF (with expecte d cosmic variance of 
~ 30% -40% at re dshifts z ~ 7-9, see Roberts onll2010blat 
Munoz et all 120101 and Section 4.2.3 of ISchenker et all 
2012al) . we anchor our constraints on the evolving UV lu- 
minosity density with precision determinatio ns of the galaxy 
lumino sity function at redshifts 4 < z < 6 by iBouwens et alj 
(l2007h . 

To utilize as much information as possible about the lu- 
minosity function (LF) constraints at z ~ 4 - 9, we perform 
Bayesian inference to generate full posterior distributions of 
the galaxy lumino sity functions at e ach redshift. To achieve 
this, we perform Schechter (1976) function parameter es- 
timation at z ~ 4 - 9 using the stepwise maximum likeli- 
hood Uy_Juminosity function constraints reported in Table 
5 of IBouwens et al. (120071 for z ~ 4 - 6) and Table 2 of 
IMcLure et al.l (120121 for z ~ 7-8) allowing all parameters 
(the luminosity function normalization (/>*, the characteristic 
galaxy magnitude M 4 , and the faint-end slope a) to vary. For 
these LF determinations, we assume Gaussian errors and a \ 2 
likelihood. For additional constraints at z ~ 7 - 8, we use the 
samples from th e full posterior distr ibutions of LF parameters 
determined by lSchenker et al.l (l2012al) . using their method for 
Bayesian inference of LF parameters discussed in their Sec- 
tion 4.2. Lastly, at redshift z ~ 9 we perform parameter esti- 
mation on given the stepwise m aximum likel i hood LF de- 
termination provided in Table 4 of lMcLure et al.l (12012) while 
keeping M+ and a fix ed at the best-fit z ~ 8 values reported by 
McLure end] dMTl . A gain we assume Gaussian errors and 
a x 2 likelihood. The limited information available at z ~ 9 
and our restricted fitting method at this redshift mean that the 
inferred allowed variation in the puy(z ~ 9) will be underes- 
timated. However, we have checked that this restriction does 
not strongly influence our results presented in Section In 
each case our maximum likelihood luminosity function pa- 
rameters ar e consistent within 1 - a o f the values originall y 
reported by IBouwens et al.l (120071) and IMcLure etail (120121) 
and are of course identical in the case of the Schenk er et al.l 



(|2012a|) LFs. 

Figure [2] shows the integrated UV luminosity density puv 
as a function of limiting magnitude for these galaxy luminos- 
ity function determinations at z ~ 5 - 8 (constraints at z ~ 4 
and z ~ 9 are also used). For each redshift, we show the max- 
imum likelihood puv (white lines) and the inner 68% varia- 
tion in the marginalized puv (blue regions). Since the lumi- 
nosity functions are steep {a < —1.7), the luminosity densi- 
ties puv increase dramatically below the characteristic mag- 
nitude M* at each redshift, but especially at z > 7. The total 
puv supplied by star-forming galaxies therefore strongly de- 
pends on the limiting magnitude adopted to limit the integral 
in Equation [5] The UDF12 campaign depth of Muv < -17 
provides puv ~ 10 26 ergs s" 1 Hz -1 Mpc~ 3 at z ~ 7, declining 
to p uv ~ 3.2 x 10 25 ergs s" 1 Hz" 1 Mpc" 3 at z ~ 8. To put 
these puv values in context, we also indicate the critical values 
of puv required to keep Qun = in Equation [T] and maintain 
reionization (grey areas and dashed line; assuming case A re- 
combinatior|3) if /esc = 0.2, log£i on = 25.2 log ergs" 1 Hz, and 
Cmi = 3. Under these reasonable assumptions, the currently 
observed galaxy population clearly is not abundant enough 
to maintain reionization at z > 7. Understanding the role of 
galaxies in the reionization process will therefore likely re- 
quire extrapolations to luminosities beyond even the UDF12 
depth, and the constraints on the LF shape achieved by the 
UDF12 observations will be important for performing these 
extrapolations reliably. The extent of the extrapolation down 
the luminosity function required for matching the reioniza- 
tion constraints depends in detail on our assumptions for the 
escape fraction / esc or the ionizing photon production rate per 
UV luminosity £j on of galaxies. We emphasize that the critical 
ionizing photon production rate n; on depends on the product 
"ion = /esc£ionPuv, and will shift up or down proportionally to 
/esc£ion at fixed puv- Similarly, the n\ on value that balances re- 
combination is proportional to the clumping factor Chii, and 
variation in the clumping factor will also shift the required 
puv up or down. 

4.1. UV Luminosity Density Likelihoods 

We wish to use the evolving UV luminosity density puv 
to constrain the availability of ionizing photons ri mn as given 
in Equation [4] Given how significantly puv increases with 
the limiting absolute magnitude Muv, we must evaluate the 
full model presented in Section [2] for chosen values of Muv- 
We select Muv = -17, Muv = -13, and Muv = -10 to provide 
a broad range of galaxy luminosities extending to extremely 
faint magnitudes. The observations probe to Muv = -17, and 
this magnitude therefore serves as a natural limit. We consider 
limits as faint as Muv = -10 as this magnitude may corre- 
spond to the minimum mass dark matter halo able to accrete 
gas from the photoheated intergalactic medium or the mini- 
mum mass dark matter halo able to retain its gas supply in the 
presence of supernova feedback. These scen arios have so far 
prove n impossible to distinguish empirically (Mu noz & Loebl 
l20Tlh . To make use of the constraints of the UV luminosity 
density on «i 0n , we need to adopt a likelihood for use with 

10 When calculating whether the production rate of ionizing photons from 
galaxies can maintain the ionization fraction of the IGM near the end of 
reionization, an assumption of case A recombination is appropriate since 
recombinations to the hydrogen ground state lar gely do not help sustai n 
the IGM ionization. For a de t ailed discussion, se e Furlanetto & Oh| 120051) ; 
Faucher-Giguere et al. (2009); Kuhlen & Faucher-Giguere (2012). We there- 
fore adopt case A recombination in Figure [2] However, we assume case B 
recombination throughout the rest of the paper. 
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FIG. 2. — Constrained ultraviolet (UV) luminosity densities puv as a function of limiting magnitude Muv and redshift t. Shown are the z, ~ 5 - 8 maximum 
likelihood values of puv vs. limiting magnitude calculated using Equation [3] (white lines), and the corresponding marginalized inner 68% credibility intervals 
for puv(^uv) (blue regions). In each panel, we indicate with a dotted line the limiting depth of the luminosity function determinations. Also shown is the 
puv required for galaxies to maintain a fully ionized universe assuming log£i on = 25.2 l og ergs" 1 Hz, / esc = 0.2, Chii = 3, and case A recombination (dashed 
lines and grey regions). We use Bayesian parameter estimation methods to d etermine th e |Schechter ( 1976) functi on parameter post erior distributions inferred 
from the stepwise maximum likelihood luminosity function (LF) data of Bouwens et al. 12007) at z. ~ 4 — 6 and McLure et al. (2012) at z. ~ 7 — 8. We also use 
the full posterior distribution sampling of the Schechter ( 1976) function parameters from the LF determination of Schenker et al. 1 2012a) to produce additional, 
independent constraints on puv at z ~ 7 — 8. At z ~ 9 where the data is limited , we si mply infer the LF normalization <j>* keeping the characteristic magnitude 
M„ and faint-end slope a fixed at the z, ~ 8 values determined by McLure et al. (2012) and expect the inferred possible variation in puv(z ~ 9) to be somewhat 
underestimated. 
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Bayesian inference on parame terized forms of A jy(z). Given 
the parameterized form of the Sch echterl (11976b function, we 
find that the marginalized posterior distributions of puv are 
skewed at z > 6, with tails extending to larger log puv val- 
ues than a Gaussian approximation would provide. To cap- 
ture this skewness, when constraining the cosmic reionization 
history in Section[5]below, we therefore use the full marginal- 
ized posterior distribution of puv provided by the integrated 
LF determinations calculated in Section|4] 

Figure [3] shows the marginalized posterior distribution of 
the UV luminosity density puv for limiting magnitudes of 
Muv = -17 (red lines), Muv = -13 (orange lines), and Muv = 
-10 (blue lin es) for our Sch e chter function fits to the z ~ 
5-6 LFs of IBouwens et all ( 120071) . the z ~ 7-8 LFs of 
ISchenker et all d2012al). and the z ~ 7-8 LFs of lMcLure etall 



(|2012|) . Although iSchenker et al.1 d2012al) and iMcLure et ail 
(2012) use the same data sets, their luminosity function deter- 
minations are based on different selection techniques. They 
therefore represent independent determinations of the high- 
redshift luminosity functions and are tre ated accordingly. We 
additionall y use constraints from z ~ 4 (IBouwens et al.ll2007h 
and z ~ 9 (Mc Lure et al.H2.Q12l) . The posterior distributions 
for puv(z ~ 9) we calculate are likely underestimated in their 
width owing to an assumption of fixed M+ and a values, but 
this assumption does not strongly influence our results. In 
what follows, these posterior distributions on puv are used as 
likelihood functions when fitting a parameterized model to the 
evolving UV luminosity density. 

5. CONSTRAINTS ON REIONIZATION 

The observational constraints on the process of cosmic 
reionization in the redshift range z > 7 are the spectral char- 
acter of high-redshift star-forming gala xies determined from 
the UDF12 program (Section |3] and iDunlop et all l2012bl) . 
the evolving luminosity densit y constraints enabled by those 
same o bser vations (Sectio n|H Ellis et al.12013 , Schenker et al. 
I2012al and lMcLure etal.ll2012l) . and the electron scattering 
optical dep th inferred from the 9-year WMAP observations of 
the CMB ilHinshaw et~aT]|2012l and briefly below). 

As we demonstrate, these constraints are in tension when 
taking the current data at face value. The Lyman contin- 
uum photon production rates per unit UV luminosity spec- 
tral density calculated using the BC03 models that are con- 
sistent with the UV spectral slopes of galaxies at z ~ 7 - 9 
are log£j on « 24.8-25.3 log ergs" 1 Hz (Section O, and 
for a reasonable escape fraction of / esc ~ 0.2 and IGM 
clumping factor of Chii, UV luminosity density of logpuv > 
26 log ergs s" 1 Hz -1 Mpc~ 3 is required to induce sig n ificant 
ionization at z > 7, but as we showed in lEllis et all d2013l) 
the observed abundance of star-forming galaxies continues 
a measured decline at high-redshift (z > 8). Further, repro- 
ducing the WMAP Thomson optical depth requires an ex- 
tended reionization process, since instantaneous reionization 
(Ghii =1) would need to occur at z = 10.3 ± 1 . 1 to repr oduce 
the measured r = 0.084 ±0.013 dHinshaw et alj 120121) and 
very likely the ionization fraction Quu < 1 at z > 7. Based 
on the UDF12 and WMAP constraints, we therefore antici- 
pate that the UV luminosity density declines to some mini- 
mum level beyond z > 7, and then persists with redshift to 
sustain a sufficient partial ionization of the IGM to satisfy the 
Thomson optical depth constraint. We need a methodology 
to quantify this process and, given its redshift dependence, to 
determine the required minimum luminosity of abundant star- 
forming galaxies to reionize the universe by z ~ 6. 



Our chosen methodology is to use a simple parameterized 
model of the evolving UV luminosity density to calculate 
the redshift-dependent ionizing photon production rate den- 
sity ni on , constrained to reproduce the UV luminosity den- 
sity values reflected by the likelihood functions shown in 
Figure [3] The evolving n; on with redshift is used to calcu- 
late the reionization history Qhu(z) by integrating Equation 
[U The corresponding Thomson optical depth is calculated 
using Equation QT| and then evaluated against the posterior 
distribution of t provided by the public 9-year WMAP Monte 
Carlo Markov Chains. At each posterior sample evaluation, 
our method requires a full reconstruction of the reionization 
history and integration of the electron scattering optical depth 
and we have incorporated the reionization calculation into 
the MultiNest Bayesian inference software (IFeroz & H obson 
120081: IFeroz et al.ll2009l) . 

5.1. A Parameterized Model for the Evolving UV Luminosity 

Density 

To infer constraints on the reionization process, we must 
adopt a flexible parameterized model for the evolving UV lu- 
minosity density. The model must account for the decline of 
puv apparent in Figures [2] and [3] without artificially extend- 
ing the trend in puv at z < 6 to redshifts z > 8. The rapid 
decline in puv between z ~ 4 and z ~ 5 suggests a trend of 
d log puv /d log z ~ -3, but the higher redshift values of puv 
flatten away from this trend, especially for faint limiting mag- 
nitudes. Beyond z ~ 10, we expect that some low-level UV 
luminosity density will be required to reproduce the Thom- 
son optical depth, to varying degrees depending on the chosen 
limiting magnitude. 

With these features in mind, we have tried a variety of pa- 
rameterized models for puv- We will present constraints using 
a three-parameter model given by 



Puv(z) = PUV, Z =4 [j£ 



+ PUV,;=7 



(?)' 



(14) 



The low-redshift amplitude of this model is anchored by 
the UV luminosity density at z ~ 4, puv.z=4 with units of 
ergs s -1 Hz -1 Mpc~ 3 . The high-redshift evolution is deter- 
mined by the normalization puv,z=7 at z ~ 7, provided in units 
of ergs s" 1 Hz -1 Mpc~ 3 , and the power-law slope 7. Using 
this model, we compute the evolving puv(z) and evaluate the 
parameter likelihoods as described immediately above. 

We have examined other models, including general broken 
(double) power-laws and low-redshift power-laws with high- 
redshift constants or fixed slope power-laws. Single power- 
law models tend to be dominated by the low-redshift decline 
in the puv and have difficulty reproducing the Thomson op- 
tical depth. Generic broken power-law models have a degen- 
eracy between the low- and high-redshift evolution, even for 
fixed redshifts about which the power laws are defined, and 
are disfavored based on their Bayesian information relative to 
less complicated models that can also reproduce the puv con- 
straints. The model in Equation [14] is therefore a good com- 
promise between sufficient generality and resulting parameter 
degeneracies. 

However, when using this model, we use a prior to limit 
7 < to prevent increasing puv beyond z ~ 9. The best cur- 
rent constraints o n the abundance of z > 8.5 galaxies is from 
lEllis et all d2013f) . where we demonstrated that the highest- 
redshift galaxies (with Muv ^ -19) continue the smooth de- 
cline in abundance found at slightly later times. Correspond- 
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FIG . 3 . — Likelihood functions of the UV luminosity density puv for different li miting magnitude s and redshifts used to constrain the reionization history. 
Shown are the z ~ 5 — 8 marginalized posterior distributions for puv determined from Schechter ( 1976) function fits to luminosity function (LF) data, as reported 
in Figureff] for limiting magnitudes of Muv < — 17 (r ed lines), Muy < —13 (orange lines), and M\jy < —10 (blue li nes). To infer the distributions at z ~ 4-9, we 
use the LF data of Bouwens et al. 1 2007) at z ~ 4-6, McLure et al. 1 2012) at z ~ 7-9, and Schenker et al. ( 2012a) at z ~ 7-8. In our posterior distributions for 
Puv(z ~ 9), we keep the characteristic magnitude M„ and faint-end-slope a values fixed at the z ~ 8 best fit values reported by McLure et al. (2012) but expect 
the width of distribution to be somewhat overly narrow. These posterior distributions on puv are used as likelihood functions when fitting a parameterized model 
to the evolving UV luminosity density. 
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ingly, the constraint on 7 amounts to a limit on introducing a 
new population of galaxies arising at z > 9. The usefulness 
of this potentiality for the reionization of the universe ha s 
been noted elsewhere (e.g.. ICenll2003t [Alvarez et al.l 120121) . 
While imposing no constraint besides 7 < 0, we typically find 
that nearly constant, low-level luminosity densities at high- 
redshift are nonetheless favored (see below). Results for mod- 
els that feature low-redshift power-law declines followed by 
low-level constant puv at high-redshift will produce similar 
quantitative results. All parameterized models we have exam- 
ined that feature a declining or constant puv and can repro- 
duce the Thomson optical depth constraint produce similar 
results, and we conclude that our choice of the exact form of 
EquationfPflis not critical. 

5.2. Reionization Constraints from Galaxies 

Figure |4] shows constraints on the reionization process cal- 
culated using the model described in Section [2] We perform 
our Bayesian inference modeling assuming Muv < -17 (close 
to the UDF12 limit at z ~ 8, maximum likelihood model 
is shown as a dashed line in all panels), Muv < -10 (an 
extremely faint limit, with the maximum likelihood model 
shown as a dotted line in all panels), and an intermediate limit 
Muv < -13 (colored regions show 68% credibility regions, 
while white lines indicate the maximum likelihood model). 
We now will discuss the UV luminosity density, stellar mass 
density, ionized filling fraction, and electron scattering optical 
depth results in turn. 

The upper left hand panel shows parameterized models of 
the UV luminosity density (as given in Equation |T4l. con- 
strained by observations of the UV luminosity density (in- 
ferred from the measured luminosity functions, see Figures 
|2] and |3]l integrated down to each limiting magnitude. The 
figure shows error bars to indicate the 68% credibility width 
of the posterior distributions of puv at redshifts z ~ 4-9 in- 
tegrated to Muv < -13, but the likelihood of each model is 
calculated using the full marginalized puv posterior distribu- 
tions appropriate for its limiting magnitude. In each case the 
Puv(z) evolution matches well the constraints provided by the 
luminosity function extrapolations, which owes to the well- 
defined progression of declining puv with redshift inferred 
from the UDF12 and earlier datasets. For reference, the max- 
imum likelihood values for the parameters of Equation [l4lfor 
Muv < -13 are log puv, : =4 = 26.50 log ergs s" 1 Hz" 1 Mpc" 3 , 
log Ajv,z=7 = 25.82 log ergs s" 1 Hz" 1 Mpc" 3 , and 7 = -0.003. 

What drives the constraints on puv for these limiting mag- 
nitudes? Some tension exists at high-redshift, as the models 
prefer a relatively flat puy(z) beyond the current reach of the 
data. In each case, the maximum likelihood models become 
flat at high-redshift and reflect the need for continued star for- 
mation at high-redshift to sustain a low-level of partial IGM 
ionization. To better answer this question, we have to calcu- 
late the full stellar mass density evolutions, reionization his- 
tories, and Thomson optical depths of each model. 

The lower left hand panel shows how the models compare 
to the stellar mas s densities extrapolated from the results by 
LStark et all (I2012T) . using the method described in Section [3] 
The error bars in the figure reflect the stellar mass densities 
extrapolated for contributions from galaxies with Muv < -13, 
but the appropriately extrapolated mass densities are used for 
each model. The error bars reflect bootstrap-calculated uncer- 
tainties that account for statistically possible variations in the 
best fit stellar mass - UV luminosity relation given by Equa- 



tion [8] and are a » 0.3 dex at z ~ 7. To balance the relative 
constraint from the stellar mass density evolution with the sin- 
gle optical depth constraint (below), we assign the likelihood 
contributions of the stellar mass density and Thomson opti- 
cal depth equal weight. The stellar mass densities calculated 
from the models evolve near the 1 - a upper limits of the ex- 
trapolated constraints, again reflecting the need for continued 
star formation to sustain a low-level of IGM ionization^ 

Assuming a ratio of Lyman continuum photon production 
rate to UV luminosity log£i on = 25.2 log ergs" 1 Hz for indi- 
vidual sourc es consistent with UV slopes measured from the 
UDF12 data dDunlop et alJ2012bl) . an ionizing photon escape 
fraction / e sc = 0.2, and an intergalactic medium clumping fac- 
tor of Chii = 3, the reionization history 2hii(z) produced by 
the evolving puv can be calculated by integrating Equation[T| 
and is shown in the upper right panel of Figure |4] We find 
that the currently observed galaxy population at magnitudes 
brighter than Muv < -17 (dashed line) can only manage to 
reionize fully the universe at late times z ~ 5, with the IGM 
50% ionized by z ~ 6 and < 5% ionized at z ~ 12. Contribu- 
tions from galaxies with Muv < -13 reionize the universe just 
after z ~ 6, in agreement with a host of additional constraints 
on the evolving ionized fraction (see Section [6] below) and 
as sug gested by some previous analyses (e.g.. lSalvaterra etafl 
1201 ll) . These galaxies can sustain a ~ 50% ionized fraction 
at z ~ 7.5, continuing to ~ 12- 15% at redshift z ~ 12. Inte- 
grating further down to Muv < -10 produces maximum like- 
lihood models that reionize the universe slightly earlier. 

The Thomson optical depths resulting from the reionization 
histories can be determined by evaluating EquationQT] To re- 
produce the 9-Year WMAP r values, most models display low 
levels of UV luminosity density that persist to high redshift. 
Even maintaining a flat puv at the maximum level allowed 
by the luminosity density constraints at redshifts z < 9, as the 
maximum likelihood models do, the optical depth just can be 
reproduced. When considering only the currently observed 
population Muv < -17, the UV luminosity density would 
need to increase toward higher z > 10 redshifts for the uni- 
verse to be fully reionized by redshift z ~ 6 and the Thomson 
optical depth to be reproduced. 

We remind the reader that while the exact results for, e.g., 
the Quii evolution of course depends on the choices for the es- 
cape fraction / esc or the ionizing photon production rate £; on . 
If /esc or £i on are lowered, the evolution of Qhii is shifted to- 
ward lower redshift. For instance, with all other assumptions 
fixed, we find that complete reionization is shifted to z ~ 5 
for / esc = 0.1 and z ~ 4.75 for log = 25.0 log ergs" 1 Hz. 
Complete reionization can also be made earlier by choosing a 
smaller Chii or delayed by choosing a larger Chii- 

Given the above results, we conclude that under our as- 
sumptions (e.g., / esc = 0.2 and Chii ~ 3) galaxies currently ob- 
served down to the limiting magnitudes of deep high-redshift 
surveys (e.g., UDF12) do not reionize the universe alone, 
and to simultaneously reproduce the puv constraints, the r 
constraint, and to reionize the universe by z ~ 6 requires 
yet fainter populations. This conclusion is a r amification of 
the UDF12 UV spectral slope constraints by iDunlop et al.1 

1 1 We note that if we ignore the electron scattering optical depth constraint, 
the other datasets favor a UV luminosity density that declines more rapidly 
at high redshift than does our maximum likelihood model calculated includ- 
ing all th e co nstraints. However, the UV luminosity density parameters in 
Equation 1 141 remain similar. Excluding the stellar mass density constraint 
has almost no effect on the maximum likelihood parameters inferred for the 
model. 
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FIG. 4. — Joint constraints on the reionization history, assuming ionizing photon contributions from galaxies with Mjjv < — 17 (maximum likelihood model 
shown as dashed line in all panels), Mijv < —10 (maximum likelihood model shown as dotted line in all panels), and M\]\ < -13 (maximum likelihood model 
shown as white line in all panels; 68% credibility regions shown as colored areas). We use the posterior distributio ns for puy with red shift shown in Figure [5] 
extrapolated stellar mass density constraints, and the posterior distribution on the electron scattering op tical depth iHinshaw et al. 2012) as likelihood functions 
to constrain the simple parameterized model for the evolving UV luminosity density given by Equation[l4]at redshifts z > 4. The constrained evolution of puv 
is shown in the upper left panel (error bars indicate the puv constraints for Muv < —13, but each model uses the appropriate constraints). From Puv(z), we 
can sim ply integrate with redshift to determine the stellar mass densities (bottom left panel; data points with error bars indicate extrapolations of the S tark et all 
12012) stellar mass densities to Muv < —13, but all models use the appropriate constraints). The models tend to exceed slightly the stellar mass densities at the 
highest-redshifts (z ~ 7), a result driven by the constraint on the election scattering optical depth. By assuming the well-motivated values of the ratio of Lyman 
continuum photon production rate to UV luminosity log£j on = 25.2 logergs -1 Hz for individual sources, an ionizing photon escape fraction / esc = 0.2, and an 
intergalactic medium clumping factor of Chii = 3, the reionization history Qhii calculated by integrating Equation[T]is shown in the upper right panel. Integrating 
the reionization history provides the electron scattering optical depth (lower right panel, 9-year WMAP constraint indicated as the grey region). 
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FIG. 5. — Reionization histories for models that include galaxies with 
M\j\ < -13 (maximum likelihood model: white line; 68% credibility re- 
gion: orange area), Mrjv < - 17 (maximum likelihood model only, dashed 
line), and Muv < — 10 (maximum likelihood model only, dotted line). We 
also show a collection of other claimed constraints on the neutral fraction 
1 _ 2hii; note how the ordinate switches from log to linear to more clearly il- 
lustrate the full range of measurements. These include (see text for details 
and references): measurements of the Lyman-a forest transmission (blue 
solid triangles), conservative upper limits on the neutral fraction from the 
fraction of dark pixels in the Lyman-a forest (purple open triangles), quasar 
near-zone measurements (open and solid magenta squares), damping wing 
absorption in a GRB (open green circle), the evolving abundance of Lyman-a 
emitter galaxies (open red pentagons), the clustering of those galaxies (filled 
dark blue diamond), and the evolving fraction of Ly man-break galaxies with 
strong Lyman-a emission lines (filled black pentagon). The grey dashed lines 
labeled "Slow kSZ" are purely illustrative, showing the slowest evolution 
currently allowed by the small-scale CMB temperature data; the endpoint of 
reionization is arbitrary in these cases. Finally, the shaded gray region shows 
the redshift during which instantaneous reionization would occur according 
to the WMAP 9-year data. 



(I2012bl) . which eliminate the possibility that increased Ly- 
man continuum emission by metal-poor p opulations could 
have produced reionization at z > 6 (e.g., Ro bertson et al.l 
2010). However, too much additional star formation beyond 
the Muv < -13 models shown in Figure |4] will begin to ex- 
ceed the stellar mass density constraints, depending on f esc . It 
is therefore interesting to know whether the Muv < -13 mod- 
els that satisfy the /?uv, p*, and t constraints also satisfy other 
external constraints on the reionization process, and we now 
turn to such an analysis. 

6. COMPARISON TO OTHER PROBES OF THE IONIZED FRACTION 

In this section we will collect constraints on the IGM neu- 
tral fraction from the literature and compare them to the evo- 
lution of this quantity in our models based on the UDF12 data. 
These constraints come from a wide variety of astrophysical 
measurements, but all are subject to substantial systematic or 
modeling uncertainties, about which we will comment below. 
We show the full set in Figure [5] along with a comparison to 
our model histories. We will first briefly discuss these con- 
straints and then how our model histories fare in comparison 
to them. 



6.1. The Lyman-a Forest 

The best kno wn Lyman-a forest constraints come from 
iFan et al.l ((2006b), who measured the effective optical depth 
evolution along lines of sight taken from SDSS (including 
both Lyman-a and higher-order transitions, where available). 
Those authors then made assumptions about the IGM temper- 
ature and the distribution of density inhomogeneities through- 
out the universe to infer the evolution of the neutral fraction 
(these assumptions are necessary for comparing the higher- 
order transitions to Lyman-a as well, because those transi- 
tions sample different parts of the IGM). The corresponding 
limits on the neutral fraction according to their model are 
shown by the filled triangles in Figure [5] 

The original data here are the transmission measurements in 
the different transitions. Transforming those into constraints 
on the ionized fraction requires a mo del that: (1) predicts 
the temperature evolution of the IGM (iHui & Haimanl 120031 : 
iTrac et al.l20 08tlFurlanett o & Ohl2009l); (2) describes the dis- 
tribut i on of gas densities in t he IGM dMiralda-Escude et al.l 
l2000tlBolton & Beckerl l2009): (3) accounts for sp atial struc- 
ture i n the averaged transmission measurements (tLidz et al.l 
2006); and (4) predicts the topology of i onized and neu- 
tral regions at the tail end of reionization (Furlanett o et al.l 
2004b; Choudh urvetaLl 120091) . These are all difficult, and 
the lFan et al.l (|2006b) constraints - based upon a simple semi- 
analytic model for the IGM structure - can be evaded in a 
number of ways. In particular, their model explicitly ignores 
the possibility that 2hii < 1, working only with the residual 
neutral gas inside a mostly ionized medium. 

Unsurprisingly, our models do not match these Lyman-a 
forest measurements very well. The crude approach of Equa- 
tion [T] fails to address effects from the detailed density struc- 
ture of the IGM on the evolution of the mean ionization frac- 
tion near the end of the reionization process. Our model is 
therefore not sufficient to model adequately the tail end of 
reionization, and only more comprehensive models that ac- 
count for both radiative transfer effects and the details of IGM 
structure will realistically model these data at z < 6. 

More useful to us is a (nearly) model-independent upper 
limit on the neutral fraction provided by sim p ly counting the 
dark pixels in the spectra (Mesinger 2010). McGreeretal. 
d201 ll) present several different sets of constraints, depend- 
ing on how one defines "dark" and whether one uses a small 
number of very deep spectra (with a clearer meaning to the 
dark pixels but more cosmic variance) versus a larger set of 
shallower spectra. Their strongest constraints are roughly 
1 - Ghii < 0.2 at z = 5.5 and 1 - g H n < 0.5 at z = 6, shown 
by the open triangles in Figure [5] Interestingly, this model- 
independent approach permits rather late reionization. 

6.2. Ionizing Background 

Another use of the Lyman-a forest is to measure the ion- 
izing background and thereby constrain the emissivity e of 
galaxies: the ionization rate T oc eA, where A is the L yman 
continuum mean free path. Bolt on & Haehneltl d2007bl) at- 
tempted such a measurement at z > 5. The method is difficult 
as it involves interpretation of sp ectra near the saturat i on limit 
of the forest. Nonetheless, the iBolton & Haehneltl (|2007b) 
analysis shows that the ionizing background falls by about a 
factor of ten from z ~ 3 to to z ~ 6 (though see McQuinn et al. 
1201 ll) . The resulting comoving emissivity is roughly the same 
as at z ~ 2, corresponding to 1 .5-3 ph otons per hydrogen atom 
over the age of the Universe (see also Haardt & Madau 2012). 
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Additional detaile d con straints were derived by 
iFaucher-Giguere et al.l (j2008), who used quasar spectra 
at 2 < z < 4 to infer a photoionization rate for intergalactic 
hydrogen. Co mbining these meas u rement s with the pre- 
vious work by iBolton & Haehnell d2007bl) and estimates 
of the int ergalactic mean free p ath of Lyman continuum 
photons dP rochaska etail 120091; ISongaila & Cowid 120101) , 
iKuhlen & Faucher-Giguerel (120121) calculated constraints 
on the comoving ionizing photon production rate. These 
inferred values of log n\ on < 51 log s" 1 Mpc~ 3 at z ~ 2-6 are 
lower than those produced by naively extrapolating typical 
models of the evolving high-redshift UV luminosity density 
that satisfy reionization constraints. To satisfy both the 
low-redshift IGM emissivity con straints and the reionization 
era co nstraints then available, IKuhlen & Faucher-G iguere 
d2012l) posited an evolving escape fraction 

/esc(z)=/0X[(l+z)/5] K . (15) 

If the power-law slope k is large enough, the IGM emis- 
sivity constraints at z < 6 can be satisfied with a low /o 
while still reionizing the universe by z ~ 6 and matching 
previous WMAP constraints on the electron scattering opti- 
cal depth. Other previous analyses have em phasized similar 
needs for an evolving escape fraction (e.g.. iFerr ara & Loeb 
l2^TllMitraet^l20L3l) . 

To study the effects of including the IGM emissivity con- 
straints by permitting an evolving / esc , we repeat our cal- 
culations in Section \5\ additionally allowing an escape frac- 
tion given by Equation Q3] with varying /o and k (and max- 
imum escape fraction / max = 1). With a limiting magni- 
tude Muv < -13, and utilizing the updated constraints from 
UDF12 and WMAP, we find maximum likelihood values 
fo = 0.054 and k = 2.4 entirely co nsistent with the results of 
Kuhlen & Faucher-Giguere (2012, see their Figure 7). When 
allowing for this evolving / esc (with / max =1), we sensi- 
bly find that the need for low-level star formation to sat- 
isfy reionization constraints is reduced (the maximum like- 
lihood model has a high-redshift luminosity density evolution 
of puv oc z~ 1,2 ) and the agreement with the stellar mass density 
constraints is improved. 

However, the impact of the evolving / esc (z) depends on the 
maximum allowed / esc . We allow only / esc to vary (not the 
product /esc^ion) owing to our UV spectral slope constraints, 
and the maximum / esc = 1 (and, correspondingly, the maxi- 
mum ionizing photon production rate per galaxy) is reached 
by z ~ 15. In this model with maximum / max = 1, the escape 
fraction at z ~ 7 - 8 (where the ionized volume filling factor 
Qhu is changing rapidly) is / esc ~ 0.17-0.22, very similar 
to our fiducial constant /esc = 0.2 adopted in Section [5] The 
rate of escaping ionizing photons per galaxy in this evolv- 
ing /esc model with maximum / max = 1 trades off against the 
high-redshift UV luminosity density evolution. The need for 
a low-level of star formation out to high redshift is somewhat 
reduced as the increasing escape fraction with redshift com- 
pensates to maintain a partial IGM ionization fraction and re- 
cover the observed WMAP t. Without placing additional con- 
straints on the end of reionization, this evolving / esc model 
with/max = 1 produces a broader redshift range of z ~ 4.5-6.5 
for the completion of the reionization process than the con- 
stant / esc model described in Section 

If we instead adopt an evolving / esc model with / max = 
0.2 we recover puv, P±, and Thomson optical depth evo- 
lutions that are similar to the constant / esc model from 



Section [5] but can also satisfy the additional low-redshift 
IGM emissivity constraints. Compared with the comov- 
ing ionizing photon production rates ri = [3.2,3.5,4.3] x 
10 50 s" 1 Mpc" 3 at redshif t s z = [4.0,4.2,5.0] inferred by 
Kuhlen & Faucher-Giguere (2012) from meas u reme nts by 
Faucher-Giguere et al.1 (120081) . iProchaska et al.l (12009b . and 
Songaila & Cowie (2010), this maximum likelihood model 
with evolving / esc recovers similar values («= [3.3,3.5,4.7] x 
10 50 s" 1 Mpc" 3 at z = [4,4.2,5]). In contrast, the constant 
/ esc model in Section [5] calculated without consideration of 
these IGM emissivity constraints would produce h ~ 8 - 13 x 
10 50 s -1 Mpc -3 at z ~ 4-5. The reduction of the escape frac- 
tion at z ~ 4-6 compared with the baseline model with con- 
stant /esc allows full ionization of the IGM to occur slightly 
later in the evolving / esc model (at redshifts as low as z ~ 5.3, 
although in the maximum likelihood model reionization com- 
pletes within Az ~ 0. 1 of the redshift suggested by the maxi- 
mum likelihood constant / esc model of Section|5]). 

With a similar r a nge of assumptions to 
IKuhlen & Faucher-Giguerel y012) for a time-dependence 
in the escape fraction, the model presented in Section [5] 
can therefore also satisfy the low-redshift IGM emissivity 
constraints without assuming an escape fraction of / esc ~ 1 at 
high red shifts. While the e volving covering fraction in galaxy 
spectra dJones et al.l 12012) provides additional observational 
support for a possible evolution in the escape fraction at 
z < 5, such empirical support for evolving / esc does not yet 
exist at z > 5. 

6.3. The Lyman-a Damping Wing 

Another use of the Lyman-a line is to measure precisely 
the shape of the red damping wing of the line: the IGM ab- 
sorption is so optically thick that the shape of this red ab- 
sorpti on wing depends upon the mean neutral fraction in the 
IGM (Miralda-Escude 1998), though the interpretation de- 
pends upon the morphology of the reion ization process, lead- 
ing to large intrinsic scatter an d biases dMcOuinn et al.ll2008t 
Mesinger & Furlanetto 2008a). Such an experiment requires 
a deep spectrum of a bright source in order to identify the 
damping wing. 

One possibility is a gamma-ray burst, which has an intrin- 
sic power-law spectrum, making it relatively easy to map the 
shape of a damping wing. The disadvantage of these sources 
is that (at lower redshifts) they almost always have damped 
Lyman-a (DLA) absorption from the host galaxy, which must 
be disentangled from any IGM signal (Chen et al. 20Q3). 

To date, the best example of such a source is GRB050904 
at z = 6.3, which received rapid followup and produced a high 
signal -to-noise spectrum (Totani et al. 2006). McOuin n et alj 
(2008) studied this spectrum in light of patchy reionization 
models. Because it has intrinsic DLA absorption, the con- 
straints are relatively weak: they disfavor a fully neutral IGM 
but allow Qhii ~ 0.5. We show this measurement with the 
open circle in Figure [5] 

Quasars provide a second possible set of sources. These are 
much easier to find but suf fer from complicated int rinsic spec- 
tra near the Lyman-a line. iSchroeder et al.l (120121) have mod- 
eled the spectra of three SDSS quasars in the range z = 6.24- 
6.42 and found th at all three are best fit if a damping wing 
is present (see also Mesinger & Haimanll2004l) . Although the 
spectra themselves cannot distinguish an IGM damping wing 
from a DLA, they argue that the latter should be sufficiently 
rare that the IGM must have Qhii 0.1 at 95% confidence. 
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We show this point as the open square in Figure This con- 
clusion is predicated on accurate modeling of the morphology 
of reionization around quasars and the distribution of strong 
IGM absorbers at the end of reionization. 

6.4. The Near Zones of Bright Quasars 

Any ionizing source that turns on inside a mostly neutral 
medium will carve out an H II region whose extent depends 
upon the total fluence of ionizing photons from the source. 
If one can measure (or guess) this fluence, the extent of the 
ionized bubble then offers a constraint on the original neutral 
fraction of the medium. This is, of course, a difficult proposi- 
tion, as the extremely large Lyman- a IGM optical depth im- 
plies that the spectrum may go dark even if the region is still 
highly-ionized: in thi s case, you find only a lowe r limit to the 
si ze of the nea r zone dBolton & Haehneltll2007al) . 

Caril li et ail (120101) examined the trends of near-zone sizes 
in SDSS quasars from z = 5.8-6.4. The sample shows a 
clear trend of decreasing size with increasing redshift (after 
compensation for varying luminosities) by about a factor of 
two over that redshift range. Under the assumption that near 
zones correspond to ionized bubbles in a (partially) neutral 
medium, the volume V oc (1 - 2hii) _1 , or, in terms of the ra- 
dius, R~^ z oc ( 1 - 2hii)- In that case, a two-fold decrease in size 
corresponds to an order of magnitud e increa se in the neutral 
fraction. However, if the Fan et al.l (1200 6b) Lyman-a forest 
measurements are correct, the neutral fraction at z ~ 5.8 is so 
small that the zones are very unlikely to be in this regime. In 
that case, the trend in sizes cannot be directly transformed into 
a constraint on the filling factor of ionized gas. We therefore 
do not show this constraint on Figure [5] as its interpretation is 
unclear. 

The rec ently discovered z = 7. 1 quasar has a ver y small 
near-zone dMortlock et al.ll201 lb . iBolton et all (1201 lb used a 
numerical simulation to analyze it in detail. They concluded 
the spectrum was consistent both with a small ionized bub- 
ble inside a region with 2hii < 0. 1 and with a highly-ionized 
medium (1 - 2hii) ~ 10 -4 — 10~ 3 if a DLA is relatively close 
to the quasar. They argue that the latter is relatively unlikely 
(occurring ~ 5% of the time in their simulations). Further 
support to the IGM hypothesis is lent by recent observations 
showing no apparent meta ls in the absorber, which would be 
unprecedented for a DLA (Si mcoe et al.H2012l) . A final possi- 
bility is that the quasar is simply very young (< 10 6 yr) and 
has not had time to carve out a large ionized region. We show 
this constraint, assuming that the absorption comes from the 
IGM, with the filled square in Figure [5] 

6.5. The Kinetic Sunyaev-Zel 'dovich Effect 

Recent small-scale temperature measurements have begun 
to constrain the contribution of patchy reionization to the ki- 
netic Sunyaev-Zel'dovich (kSZ) effect, which is generated by 
CMB photons scattering off coherent large-scale velocities in 
the IGM. These scatterings typically cancel (because any red- 
shift gained from scattering off of gas falling into a poten- 
tial well is canceled by scattering off gas on the other side), 
but during reionization mo dulation by the ionizat i on field can 
preve nt such cancellation (Gruzin ov & Hul [T998; Knox et al. 
|1998|). There is also a contribu tion from nonlinear evolution, 
see lOstriker & Vishniac] (fl986b . Both the Atacama Cosmol- 
ogy Telescope and the So uth Pole Telescope have placed up- 
per lim its on this signal (Dun klev et al.ll201lt iReichardt et al.l 
12011) . 



Because the patchiness and inhomogeneity of the process 
induce this signal, the kSZ signal grows so long as this patchy 
contribution persists. As a result, it essenti ally constrains the 
duration of reionization. Using SPT data, IZahn et ail (120121) 
claim a limit of Az < 7.2 (where Az is the redshift differ- 
ence between (2hii = 0.2 and £>hii = 0.99) at 95% confidence. 
The primary limiting factor here is a potential correlation 
between the thermal Sunyaev-Zel'dovich effect and the cos- 
mic infrared background, which pollutes the kSZ signal. If 
that possib ility is ignored, the lim it on the duration falls to 
Az < 4.4. Mesinge r et al.l (120 121) found similar limits from 
ACT data. They showed that the limit without a correlation 
is difficult to reconcile with the usual models of reionization. 
Intensive efforts are now underway to measure the correlation. 

We show two examples of the slowest possible evolution 
allowed by these constraints (allowing for a correlation) with 
the straight dashed lines in Figure [5] Note that the timing 
of these curves is entirely arbitrary (as is the shape: there is 
absolutely no reason to expect a linear correlation between 
redshift and neutral fraction). They simply offer a rough guide 
to the maximum duration over which substantial patchiness 
can persist. 

6.6. Lyman-a Lines in Galaxies 

The final class of probes to be considered here relies on 
Lyman-a emission lines from galaxies (including some of 
those cataloged in the HUDF09 and UDF12 campaigns). As 
these line photons propagate from their source galaxy to the 
observer, they pass through the IGM and can suffer absorption 
if the medium has a substantial neutral fraction. This absorp- 
tion is generally due to the red damping wing, as the sources 
most likely lie inside of ionized bubbles, so the line photons 
have typically redshifted out of resonance by the time they 
reach neutral gas. 

One set of constraints come from direct (narrowband) 
surveys for Lyman-a emitting galaxies. As the IGM be- 
comes more neutral, these emission lines should suffer 
more and more extinction, and so we should see a drop 
in the number density of objects se l ected in this manner 
(ISantosI 12001 iFurlanetto et al.ll2004al l2006t iMcOuinn efafl 
2007; Mesinger & Furlanetto 200 83). Such surveys have a 
long history (e.g. [Hu et al.ll2002t iMalhotra & Rhoadsl 120041 
ISantos et al]|2004t iKashikawa et alj|2006l) . To date, the most 
comprehens i ve sur veys have come from the Subaru telescope. 
lOuchi et"ai] (120101) examined z = 6.6 Lyman-a emitters and 
found evidence for only a slight decline in the Lyma n-a trans- 
missio n. They estimate that Qnn > 0.6 at that time. lOta et alJ 
(2008) examined the z = 7 window. Very few sources were 
detected, which may indicate a rapid decline in the popula- 
tion. They estimate that Quu ~ 0.32-0.64. We show these 
two constraints with the open pentagons in Figure [5] 

The primary difficulty with measurements of the evolu- 
tion of the Lyman-a emitter number density is that the over- 
all galaxy population is also evolving, and it can be diffi- 
cult to determine if evolution in the emitter number counts 
is due to changes in the IGM properties or in the galaxy 
population. An alternate approach is therefore to select a 
galaxy sample independently of the Lyman-a line (or at least 
as independently as possible) and determine how the frac- 
tion of these gal axies with strong Lyman-a lines evolves 
dStark et al.ll2010t) . Such a sample is available through broad- 
band Lyman break searches with HST and Suba ru. Several 
groups have performed these searches to z ~ 8 ( Fontan a et al. 
I2010t [Pentericc iet~aT1l201 U ISchenker et"ai1l2012bllOno et al. 
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12012b!) and found that, although the fraction of Lyman-a 
emitters increases slightly for samples of similar UV lumi- 
nosity fro m z ~ 3-6, there is a marked decline beyond z ~ 6.5 
(see also iTreu et al.l l20T2h . This decline could still be due 
to evolution ary processes within t he population (e.g., dust 
content, see iDayal & Ferraral l2012). but both the rapidity of 
the possible evolution and its reversal from trends now well- 
established at lower redshift make this unlikely. Assuming 
that this decline can be attributed to th e increasing neutral- 
ity of the IGM, it requires Qhit < 0-5 (jM cOuinn et al] |2007t 
Mesinger & Furlanettol20 08b: Diik stra et al.l201 ll) . We show 
this constraint with the filled pentagon in Figure [5] 

There is one more signature of reionization in the Lyman- a 
lines of galaxies. A partially neutral IGM does not extinguish 
these lines uniformly: galaxies inside of very large ionized 
bubbles suffer little absorption, while isolated galaxies disap- 
pear even when Quu is large. This manifests as a change in 
the apparent clustering of the galaxies, which is attractive be- 
cause such a strong change is difficult to mimic with baryonic 
processes within and around galaxi es (Furla netto et al.ll2006t 
iMcOuinn etail 120071: iMesinger & Furlanetto! I2008bl) . Clus- 
tering is a much more difficult measurement than the num- 
ber density, requiring a large number of sources. It is not 
yet possible with the Lyman-a line sources at z ~ 7, but at 
z ~ 6.6 th ere is no evidence for an anomalous in crease in 
clustering (McOuin net al.ll2007t lOuchiet al]|2010t) . indicat- 
ing Quu > 0.5 at that time. We show this constraint with the 
filled diamond in Figure [5] 

6.7. Comparison to Our Models 

Figure|5]also shows the reionization history in our preferred 
model, with the associated confidence intervals, that extrapo- 
lates the observed luminosity function to Muv = -13. Amaz- 
ingly, this straightforward model obeys all the cons traints we 
have lis ted in this section, with the exception of the iFan et al] 
(2006b) Lyman-a forest measurements. However, we remind 
the reader that our crude reionization model fails at very small 
neutral fractions because it does not properly account for the 
high gas clumping in dense systems near galaxies, so we do 
not regard this apparent disagreement as worrisome to any de- 
gree. 

It is worth considering the behavior of this model in some 
detail. A brief examination of the data points reveals that the 
most interesting limits come from the Lyman-a lines inside 
of galaxies (the filled pentagon, from spectroscopic followup 
of Lyman-break galaxies, and the open pentagons, from direct 
narrowband abundance at z ~ 6.6 and 7). These require rela- 
tively high neutral fractions at z ~ 7 in order for the IGM to 
affect the observed abundance significantly; the model shown 
here just barely satisfies the constraints. In a model in which 
we integrate the luminosity function to fainter magnitudes 
(e.g., Myv < -10, dotted line), the ionized fraction is some- 
what higher at this time and may prevent the lines from being 
significantly extinguished, while a model with the minimum 
luminosity closer to the observed limit (e.g., Muv < -17, 
dashed line) reionizes the universe so late that the emitter pop- 
ulation at z ~ 6.6 should have been measurably reduced. 

Clearly, improvements in the measured abundance of strong 
line emitters at this epoch (perhaps with new multi-object in- 
frared spectrographs, like LUCI on LBT or MOSFIRE on 
Keck, or with new widefield survey cameras, like Hyper- 
SuprimeCam) will be very important for distinguishing viable 
reionization histories. Equally important will be improve- 
ments in the modeling of the decline in the Lyman-a line 



emitters, as a number of factors both outside (the morphology 
of reionization, resonant absorption near the source galaxies, 
and absorption from dense IGM structures) and inside galax- 
ies (dust absorption, emission geometry, and winds) all af- 
fect the detailed interpretation of the raw measurements (e.g ., 
ISantosll2004t lDijkstraetaD20l U iBolton & Haehneltl feO 1 2) . 

In order to have a relatively high neutral fraction at z ~ 7, a 
plausible reionization history cannot complete the process by 
z ~ 6.4 (to which the most distant Lyman-a forest spectrum 
currently extends), unless either: (1) the Lyman-continuum 
luminosity of galaxies, per unit 1500 A luminosity, evolves 
rapidly over that interval, (2) the escape fraction / esc increases 
rapidly toward lower redshifts, or (3) the abundance of galax- 
ies below Muv ~ -17 evolves rapidly. Otherwise, the UDF12 
luminosity functions have sufficient precision to fix the shape 
of the Qun(z) curve over this interval (indeed, to z ~ 8), and 
as shown in Figure [5] the slope is relatively shallow. Indeed, 
in our best-fit model the universe still has 1.0— Quu ~ 0.1 at 
z ~ 6. This suggests intensifying searches for the last neutral 
regions in the IGM at this (relatively) accessible epoch. 

Although reionization ends rather late in this best-fit model, 
it still satisfies the WMAP optical depth constraint (at I — a) 
thanks to a long tail of low-lev el star formation o ut to very 
high redshifts, as suggested by lEllis et al. (2013), although 
the agreement is much easier if the WMAP value turns out 
to be somewhat high. This implies that the kinetic Sunyaev- 
Zel'dovich signal should have a reasonably large amount of 
power. Formall y, our best-fit mode l has Az ~ 5, according to 
the definition of Zah n et al.l (120121) . which is within the range 
that can be constrained if the thermal Sunyaev-Zel'dovich 
contami nation can be sorted o ut. We also note that the WMAP 
results ( Hinsha vv et al.l I20L2T) indicate that when combining 
the broadest array of available data sets, the best fit electron 
scattering optical depth lowers by ~ 5% tor « 0.08. 

In summary, we have shown that the UDF12 measurements 
allow a model of star formation during the cosmic dawn 
that satisfies all available constraints on the galaxy popula- 
tions and IGM ionization history, with reasonable extrapola- 
tion to fainter systems and no assumptions about high-redshift 
evolution in the parameters of star formation or UV photon 
production and escape. Of course, this is not to say such 
evolution cannot oc cur ("see lKuhien & F aucher-Giguer d2012l 
lAlvarez et al.l2012l and our Section l6~51 for examples in which 
it does), but it does not appear to be essential. 

7. SUMMARY 

The 2012 Hubble Ultra Deep Field (UDF12) campaign, 
a 128-orbit Hubble Space Telescope (HST) program (GO 
12498, PI: R. Ellis, a s described in lEllis et alJl20T3l and 
iKoekemoer et alj 120131) . has acquired the deepest infrared 
WFC3/IR images ever taken with HST. These observations 
have enabled the first identifi cation of galaxie s at 8.5 < z < 
12 in the Ultra Deep Field (lEllis et al.ll2013l) . newly accu- 
rate luminosity function determinations at r edshifts z ~ 7 - 8 
dSchenker et al.ll2012at iMcLure et al.l 120121) . robust determi- 
nations of the ultravio let spectral slopes of galaxies at z ~ 7-8 
(Dunl opet al.ll2012bl) . and the first esti mates of the luminos- 
ity function and spec tral slopes at z ~ 9 (IMcLure et alJl2012t 
iDunlop et al.ll201 2b). Synthesizing these constraints on high- 
redshift galaxy populations with the recent 9-year Wilkin- 
son Microwave Anisotropy Probe (WMAP) constraints on the 
electron scatter ing optical depth g enerated by ionized cos- 
mic hydrogen (Hinshaw et al. 2012) and previous determina- 
tions of the galaxy luminosity function at redshifts 4 < z < 6 
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dBouwens et al.l 120071) . we infer constraints on the reioniza- 
tion history of the universe. 

First, we use the U V spectral slope = —2 of high-redshift 
galaxies measured bv iDunlop et al.l (12012b!) to constrain the 
available lBruzual & Charlotl (120031) stellar population models 
consistent with the spectral character of galaxies at z ~ 7-9. 
These models motivate the adoption of a Lyman continuum 
photon production rate per unit UV luminosity spectral den- 
sity log£i on = 25.2 log ergs" 1 Hz; the data does not favor 
a luminosity-dependent variation in the efficiency of Lyman 
continuum photon production. With this value of £i on for high- 
redshift galaxies, and under reasonable assumptions for the 
Lyman continuum photon escape fraction from galaxies and 
the clumping factor of intergalactic gas (as motivated by cos- 
mological simulations), we find that the currently observed 
galaxy population accessible to the limiting depth of UDF12 
(Muv < -17 to z ~ 8) cannot simultaneously reionize the uni- 
verse by z ~ 6 and reproduce the Thomson optical depth r 
unless the abundance of star-forming galaxies or the ionizing 
photon escape fraction increases beyond redshift z ~ 12 from 
what is currently observed. 

If we utilize constraints on the evolving galaxy luminosity 
function at redshifts 4 < z < 9 to extrapolate down in luminos- 
ity, we find that the tension between the declining abundance 
of star-forming galaxies and their stellar mass density with 
redshift, the observed requirement to reionize the universe by 
z ~ 6, and reproducing the large electron scattering optical 
depth r « 0.084 is largely relieved if the galaxy population 
continues down to Muv < -13 and the epoch of galaxy forma- 
tion continues toz^l2-15. Given the first identification of a 



z ~ 12 candidate galaxy by the UDF12 program, the prospect 
for high-redshift galaxies to reionize the universe is positive 
provided that the epoch of galaxy formation extends to z > 12. 
Further observations by HST (e.g., the "Frontier Fields") and, 
ultimately, the James Webb Space Telescope will be required 
to answer these questions more definitively. 
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